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ABSTRACT 

We study the X-ray emission in a sample of galaxy clusters using the BeppoSAX PDS instrument in the 
20 - 80 keV energy band. We estimate the non-thermal hard X-ray cluster emission (HXR) by modeling 
the thermal contribution from the cluster gas and the non-thermal contamination from the unobscured 
AGN in the clusters. We also evaluate the systematic uncertainties due to the background fluctuations. 
Assuming negligible contamination from the obscured AGN, the resulting non-thermal component is 
detected at a 2a level in ^50% of the non-significantly AGN-contaminated clusters: A2142, A2199, 
A2256, A3376, Coma, Ophiuchus and Virgo. The data are consistent with a scenario whereby relaxed 
clusters have no hard X-ray component of non-thermal origin, whereas merger clusters do, with a 20 - 
80 keV luminosity of ^ iQ-i3-4:4 ]^-2 ^-i^ r^^^ co-added spectrum of the above clusters indicates a 

power-law spectrum for the HXR with a photon index of 2.8j;Q'4 in the 12 - 115 keV band, and we find 
indication that it has extended distribution. These indications argue against significant contamination 
from obscured AGN, which have harder spectra and centrally concentrated distribution. These results 
are supportive of the assumption of the merger shock acceleration of electrons in clusters, which has 
been proposed as a possible origin of the non-thermal hard X-ray emission models. Assuming that 
the Cosmic Microwave Background photons experience Inverse Compton scattering from the merger- 
accelerated relativistic electrons, and thus produce the observed HXR, the measured hard X-ray slope 
corresponds to a differential momentum spectra of the relativistic electrons with a slope of /i = 3.8 - 5.0. 
In presence of cluster magnetic fields this relativistic electron population produces synchrotron emission 
with a spectral index of 1.4 - 2.1, consistent with radio halo observations of merger clusters. Thus both 
hard X-ray and radio observations of merger clusters are consistent with the Inverse Compton model. 
The observed slope of HXR is also consistent with that predicted by the non-thermal bremsstrahlung, 
which thus cannot be ruled by the fit to the current data, even though this model requires an extreme, 
untenable cluster energetics. Assuming centrally concentrated distribution of HXR, the data requires a 
harder slope for the HXR spectrum, which is consistent with secondary electron models, but this model 
yields a worse fit to the PDS data and thus seems to be disfavored over the primary electron Inverse 
Compton model. 

Subject headings: galaxies: clusters - X-rays: galaxies - radiation mechanisms: non-thermal 



1. INTRODUCTION 

Non-thermal hard X-ray (HXR) emission has recently 
been observed in several clusters and groups of galaxies 
with the MECS and PDS instruments onboard BeppoSAX 
and with GIS instrument onboard ASCA. In the case of 
Coma (Fusco-Femiano et al. 1999), A2256 (Fusco-Femiano 
et al. 2000), HCG62 (Fukazawa et al. 2001) and A754 
(Fusco-Femiano et al. 2003) detection of excess emission 
above the contribution from the hot ICM is statistically 
significant, while marginal evidence is provided for A3667 
(Fusco-Femiano et al. 2001) and A2199 (Kaastra et al. 
1999). 

Most models of the HXR emission involve a population 
of electrons accelerated in the cluster medium. A natu- 
ral source of acceleration in clusters is provided by merger 
shocks. In a strong cluster merger event, the electrons 
are accelerated to relativistic speeds (e.g. Bell, 1978a, b; 
Fujita & Sarazin 2001; Takizawa & Naito 2000). The 
Inverse Compton scattering of Cosmic Microwave Back- 
ground photons from the relativistic electrons in clusters 



then can produce a non-thermal tail which exceeds the 
thermal bremsstrahlung emission at energies above 20 keV 
(e.g. Sarazin 1999, Blasi & Colafrancesco 1999). This 
model has been proposed as the simplest possible expla- 
nation of the HXR properties of galaxy clusters. If the ac- 
celeration is provided by a less energetic merger or turbu- 
lence (e.g. Ensslin et al. 1999), or if there is a high-energy 
cutoff in the electron velocity distribution, the resulting 
electron population is effectively transrelativistic. In this 
case, the dominating mechanism in producing hard X-rays 
has been proposed to be non-thermal bremsstrahlung (e.g. 
Sarazin & Kempner 2000). However, this solution faces 
some crucial problems mainly concerning the large energy 
injection required by such a mechanism and the resulting 
large heating expected (e.g., Petrosian 2001). 

In the secondary electron population models the merger 
shocks and galaxy activity accelerate and inject large 
quantities of relativistic protons into the cluster atmo- 
sphere. Most of the relativistic protons can be confined 
and accumulated in the cluster medium for very long 
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times, comparable with the cluster age ~ , and can 
then produce secondary electrons via proton-proton col- 
lisions (e.g., Colafrancesco and Blasi 1998). The energy 
losses are balanced by the continuous refilling of the new 
electrons produced in situ (i.e. continuously in time and 
everywhere in space) . Thus the resulting HXR spectrum in 
the secondary models reflects the electron spectrum right 
after the acceleration event, while in the primary models 
the more energetic electrons loose energy rapidly and the 
HXR spectrum steepens accordingly with time. 

In the present work, we expand the database of cluster 
hard X-ray emission by studying a sample of clusters ob- 
served with the BeppoSAX PDS. We model the thermal 
and AGN contributions in the sample in order to obtain 
estimates for the non-thermal component. We propagate 
the modeling uncertainties, as well as the background fluc- 
tuation uncertainties in order to obtain reliable (and some- 
what conservative) estimates for the non-thermal compo- 
nent. We furthermore study the co-added non-thermal 
hard X-ray spectrum of the sample, in order to investigate 
the origin of this emission. 

We consider uncertainties and significances at 1 cr level, 
and use H = 50 x hso km s~^ Mpc~^, unless stated oth- 
erwise. We define HXR in this paper as 20 - 80 keV PDS 
net count rate, after removing the sky background, cluster 
thermal component and AGN contamination. 

2. PDS ANALYSIS 

The sample consists of all publicly available clusters (as 
of June 2001) observed with BeppoSAX whose tempera- 
tures are constrained within ^10% . 

2.1. Data processing 

The observations were processed using SAXDAS 2.2.1. 
Extreme care was taken during the processing of the PDS 
data. We removed spikes which are caused by charged par- 
ticles hitting only one of the collimators using the method 
described by Guainazzi et al. (1997). The effect of the 
spike-removal was negligible compared to the statistical 
uncertainties. 

For A1795, A2163, A2256, A3667 and Coma the obser- 
vations were divided into several exposures, and A3627 has 
several available pointings. The datasets were processed 
separately, and the resulting spectra were co-added. 

In order to improve the S/N we restricted the PDS anal- 
ysis to 20-80 keV band and binned the data to contain a 
single bin covering this band. 

2.2. Background subtraction 

The total background in 20-80 keV band is ~10 c/s 
(Frontera et al. 1997b). Because the clusters in our sample 
have lower count rates than the background, we addressed 
carefully the uncertainties involved in the background sub- 
traction. In the standard observing mode, the PDS system 
of two identical collimators is rocked back and forth after 
each 96s, keeping one collimator at all times pointed at 
the X-ray target, to allow the simultaneous monitoring of 
source and background (Frontera et al. 1997a). In the 
standard data processing, both offset-positions (3.5° away 
from the source) are used for the background subtraction. 
The allowed upper limit for the background modulation 
with the rocking collimator offset angle is only 2% (Fron- 



tera et al. 1997a). Possible variations in the cosmic ray in- 
duced internal background due the changing environment 
are addressed by the standard method of subtracting the 
simultaneous background spectrum obtained from regions 
close to the source. The hard X-ray sky background is 
composed mostly of discrete sources (absorbed AGNs) and 
is not uniform (e.g. Comastri et al. 1995). Fluctuations 
in the number and flux of the sources could give rise to a 
significant uncertainty in the background subtraction. Re- 
lated to this is the effect of the presence of weak sources in 
only one of the offset-positions which is used for the back- 
ground subtraction (the same effect as described above, 
except that the strength of the source is such that it is 
discernible in one pointing). In order to estimate the ef- 
fect of random faint sources (e.g. weak AGNs) in the FOV 
in the offset positions, we determined the count rates for 
the cluster using only one of the two (either the negative or 
positive direction). A clear difference between the spectra 
could indicate the presence of source (s) in the offset po- 
sition. Note that when using only one offset position for 
the background the exposure time of the background is 
effectively halved, which results in a larger uncertainty in 
the source fluxes. 

Initially, in order to test the robustness of the stan- 
dard method in the case of no background fluctuations, 
we used a small sample of 12 pointings (Polaris, Lockmann 
hole, secondary pointings) in which no sources (other than 
"background") are thought to be present. From this we 
concluded that the mean 20-80 keV flux is very close to 
zero (-(0.4±3.7)10~^ c s~^) when using the standard back- 
ground subtraction, indicating that the systematic uncer- 
tainty of the standard background subtraction is negligible 
compared to statistical errors (Table However, when 
using either the positive or negative offset directions only 
we noticed a systematic difference of ~Q.06 counts s^^ be- 
tween these two with the positive pointings giving higher 
background subtracted source count rates. 

In order to study the above difference due to different 
offset pointings in more detail we selected a sample of 164 
pointings for which the (non-spike filtered) exposure times 
were larger than 20 ks (in order to have similar exposures 
to the cluster sample) and for which the count rate ob- 
tained with the standard processing (including spike fil- 
tering) is between -0.1 and 0.1 counts s~^. These crite- 
ria were chosen in order to select either blank fields, or 
faint sources. The mean exposure time ranges between 
19.1 ks and 78.3 ks with a mean of 32.2 ks. The spa- 
tial distribution of the analyzed pointings is rather uni- 
form. We find (see Fig. that the mean count rate as 
obtained through the standard analysis is 0.027 ± 0.051 
counts s^^, which suggests a non-significant detection in 
the whole sample. The difference between the count rates 
obtained with either the negative or positive offset direc- 
tion amounts to 0.058 counts s~^ and the distributions 
are symmetric about the mean obtained with the standard 
analysis. We therefore use a correction of 0.029 counts s~^ 
(with the proper sign) whenever we use the count rates ob- 
tained using only one offset position. A possible reason for 
this effect could either be the effect of radiation entering 
the collimators from the side, screening of the instruments 
by the satellite, or the fact that the detector is looking at 
more/less radioactive parts of the satellite. 

Furthermore we quantified the effect of the background 
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fluctuations. To do this, we used the above sample of 164 
pointings to compute the differences in background sub- 
tracted count rates between positive and negative back- 
ground pointing directions. We divided them into 4 ex- 
posure time bins (23, 35, 44, 60 ks) and determined the 
widths (tr) of the obtained distributions. The widths de- 
crease with increasing exposure time, as expected since 
the widths of the distributions can be described by cr^ = 
"■stat + '''flue I where CTstat is the statistical uncertainty (de- 
pendent on the total number of counts in the background 
or the exposure time), and cr^^^. is due to the real back- 
ground variations. We then assumed that crgtat is propor- 
tional to \fi (where t is the exposure time) and fitted this 
function, obtaining crfluc = 0.027 cts s~^. In the standard 
background subtraction the uncertainty introduced by the 
background fluctuations should be lower by a factor \/2, 
i.e., 0.019 cts s^^, since 2 background fields are used. In 
the following analysis we use this value as a systematic er- 
ror in the background subtracted source count rates when 
using the standard method and the above value 0.027 cts 
when only one background pointing direction (positive 
or negative) is used. 

In Fig. 121 we plot the difference between the count rates 
of the two different background pointings for all cluster 
exposures. The data give an average difference of 0.068 
c s""'^, very close to the value found in the blank fields 
above. Eight exposures deviate from the mean at 90% 
confidence level, while random fluctuations would predict 
only 4, thus our sample is likely contaminated by point 
sources in background regions. Thus, we reject the >90% 
deviant pointings for clusters A85, A1750, A2142, A2390, 
A3562, A3571, and RXJ1347.5-1145, and correct the re- 
sulting data values by the systematic shift ±0.029 c s^^ 
found above using the blank fields (for A2163 only a frac- 
tion of the total exposure is affected and thus we make 
no correction for it). This removes most of the negative 
count rate detections, with a notable exception of A3571. 
The signal of A3571 in 20-80 keV band remains negative 
regardless of whether we use one offset pointing or two, 
indicating that possibly both background offsets are con- 
taminated by AGNs. We chose to use the positive offset for 
background subtraction for A3571, because this gives a net 
count rate closest to zero and thus minimizes the oversub- 
traction due to AGNs. For the rest we use the standard 
method of using both positive and negative background 
pointings. We add the systematic errors due to the back- 
ground fiuctuation found above in quadrature to the sta- 
tistical uncertainties of each cluster. The obtained count 
rates are listed in Table ^ Using the total 20 - 80 keV 
band emission, we achieve 3cr detections in the direction 
of the following 10 clusters: A2142, A2256, A3376, A3627, 
A3667, Coma, Cygnus A, Ophiuchus, Perseus and Virgo. 

2.3. Vignetting 

The vignetting of PDS is assumed not to vary with pho- 
ton energy (Frontera et al. 1997). It is modeled with 
a linear function of off-axis angle, reaching zero at 1.3°. 
When predicting the thermal contribution at 20 - 80 keV 
energies, we use the actual PDS data (see Section|3} at 12 
- 20 keV to normalize the model and thus the vignetting 
is taken care of without further work. When predicting 
the AGN contribution, we correct it by multiplying by the 
vignetting factor at the off-axis of the AGN. When de- 



riving the non-thermal hard X-ray luminosity for a given 
cluster, we assume that the HXR is distributed like the 
intracluster gas and thus we use the vignetting function 
and the /3 model to obtain the vignetting correction to the 
luminosity obtained with the on-axis response. The effect 
of the vignetting correction to the luminosity is small, at 
the level of at most 20% for the closest clusters, and thus 
the assumption on the spatial distribution of HXR is not 
important. 

3. THERMAL MODELS 

We model the thermal component with XSPEC model 
wabs X mekal. The BeppoSAX MEGS study of most of 
these clusters reports cluster average temperatures (de- 
Grandi et al. 2001). Since in that work the central regions, 
affected by the presence of cooler gas (the "cooling fiow" 
scenario) are excised from the estimates while they are 
included in PDS data, we preferred to use the published 
ASCA single temperature fits (Markevitch et al. 1998) for 
the "cooling flow" clusters, and BeppoSAX values for the 
"non-cooling flow" clusters. For the distant (z>0.1) clus- 
ters only BeppoSAX results are available and thus we use 
these regardless of the presence of any cooler gas. 

MEGS and ASGA do not cover the fuU FOV of PDS. 
Thus, due to the radially decreasing temperature pro- 
flies consistently observed with BeppoSAX (deGrandi & 
Molendi 2002) and ASGA (Markevitch et al. 1998), the 
temperatures obtained by these instruments are high com- 
pared to the global ones. For this reason we compare our 
adopted temperatures with results obtained for a subsam- 
ple, with Ginga, because it covers the full FOV of PDS. 
Ginga temperatures for A496, A 179 5, A2142 and A2199 
(White et al. 1994) are systematically, and in most cases 
signiflcantly, below the ASGA values, consistently with 
the radially decreasing temperature proflles. However, the 
thermal model predictions in PDS 20 - 80 keV band us- 
ing either values are consistent within the thermal model 
normalization uncertainties and thus in general the radial 
temperature decrease in clusters does not significantly af- 
fect the conclusions of this work. 

For the nearby clusters Goma, Ophiuchus, Perseus and 
Virgo, the fraction of the cluster covered by MEGS and 
ASGA is small and thus the effect of decreasing tempera- 
ture profiles may be significant for these. For example, the 
- 1° FOV Ginga temperatures of Goma (8.21±0.16 keV, 
Hughes et al. 1993) and Perseus (6.33tQ^8, Allen et al, 
1992) are smaller than the corresponding MEGS 0-20' val- 
ues (9.20±0.13 keV deGrandi et al. 2002 and 6.68±0.08, 
respectively). The model prediction in 20 - 80 keV band 
of PDS with Ginga parameters for Goma is significantly 
smaller than with BeppoSAX values. Thus, for the nearby 
Goma and Perseus clusters we adopt the Ginga results. For 
the Virgo cluster Ginga data are not available and we use 
the MEGS to estimate the temperature (see Section 0J. 
For Ophiuchus there are ^ 1° FOV Tenma results (11.6 
keV; Okumura et al. 1998) but these are at odds with 
the 0'-8' MEGS value of 10.9±0.3 keV. Using the Tenma 
value, and normalizing the model to PDS 12 - 20 keV (see 
below) , the model prediction is significantly above the ob- 
served emission, indicative of overestimation of the tem- 
perature. For a hot and nearby cluster like Ophiuchus, 
the PDS data are of sufficiently good quality for the pur- 
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pose of spectral fitting. In the Opliiuclius field there are 
no contaminating AGN (see Section 0J and thus we can 
assume that the low energy band of PDS (12-35 keV) is 
dominated by the thermal emission of the whole cluster. 
Fitting this band with mekal, keeping the metal abun- 
dance fixed to 0.3 Solar, we obtained a temperature of 
9.1±0.6 keV. MEGS and PDS values are consistent with 
the decreasing temperature profile (and inconsistent with 
Tenma values), and thus we adopt our PDS results for the 
thermal model of Ophiuchus. 

In order to normalize the above models to the larger (ra- 
dius of 1.3°) FOV of PDS, we fitted the thermal models 
to the 12 - 20 keV band PDS data with the normalization 
as the only free parameter. To check the robustness of the 
fit, we predicted the normalization using /3 models to com- 
pute the increment of the model normalization between 
the region where the cluster model has been normalized, 
usually by ROSAT PSPG (Ebeling et at. 1996), and the 
PDS FOV. The fitted and predicted normalizations dif- 
fer by more than 50% for the faint clusters (A348, A1750, 
A2390, PKS0745-191, RXJ0152.7-135.7 , RXJ1347.5-1145 
and Zw3146), probably due to large statistical uncertain- 
ties in the 12-20 keV data. Thus we reject these clusters 
from further analysis. For the rest these two methods give 
values that differ by less than 40%, which can be explained 
by the uncertainties involved in the radial extrapolation 
and the cross-calibration uncertainty between PDS and 
PSPG. We prefer to use the fitted values, because they 
should be devoid of these uncertainties. 

In Virgo, Perseus and Gygnus A the 2-10 keV band data 
are strongly contaminated by AGNs. Thus, when normal- 
izing the thermal model, the AGN contribution must be 
taken into account. We will describe this in detail in Sec- 
tion 4. 

As a further check of the bright AGNs in background 
fields, we compared the obtained thermal model normal- 
izations in 12 - 15 and the 15 - 20 keV energy bands. The 
appearance of an AGN in only one background field would 
result in different normalization using either positive or 
negative pointing for background subtraction. Also, if the 
spectrum of an AGN is not identical to that of the clus- 
ter in the 10-20 keV band, its appearance would result 
in different normalization using either the 12-15 or 15-20 
keV band. We found that in the 15-20 keV band in all 
clusters, except A3571, both offsets give consistent values 
for the normalization of the thermal model. In the 12-15 
keV band clusters A3571 and A2142 give inconsistent val- 
ues between the two offsets. Therefore, for A2142 we use 
only the 15-20 keV band for the normalization, and for 
A3571 we use the predicted normalization. 

The obtained model predictions are listed in Table 1. 
The reported uncertainties of the thermal models include 
only the statistical uncertainty due to the PDS data in 
12-20 keV band. In most cases this is negligible com- 
pared to the PDS 20-80 keV count rate uncertainties. We 
previously discovered that for a subsample of clusters pos- 
sible hot IGM temperature variations result in negligible 
variation of the thermal model prediction in the 20 - 80 
keV band. Assuming that this holds for the whole sam- 
ple, propagating the model normalization uncertainty only 
is adequate to estimate the uncertainties in the thermal 
model. 



4. AGN 

In the large field of view of PDS the 20-80 keV band 
emission may be contaminated by AGN and QSO ran- 
domly projected in the line of sight, or AGN belonging to 
the cluster under study. In optical surveys it was found 
that most of the nearby AGN (^ 80%) are optically faint 
Scyfert 2 galaxies (e.g. Maiolino & Rieke 1995). In the uni- 
fied AGN scheme the optically bright and identifiable AGN 
or Seyfert 1 (Syl) are the ones observed face-on with no 
obscuration by the torus. Most of the lines of sight to the 
AGN nucleus are intersected by the absorbing torus and 
thus most of the AGN are obscured (NH = 10^2-25 a^oms 
cm~^, Risaliti et al. 1999) and optically faint Seyfert 2 
(Sy2). Recent deep X-ray observations of blank fields (e.g. 
Hasinger et al. 2001) have consistently discovered a pop- 
ulation of absorbed point sources which outnumbers the 
Syl by a factor of ~ 4. In addition, a population syn- 
thesis modeling of Cosmic X-ray Background (Gilli et al., 
1999) indicates that 80% of the AGN need to be obscured 
to produce the GXB spectrum, which is harder than the 
spectrum of unobscured AGN. 

The local background has been subtracted from the PDS 
data and thus the effect of random AGN and QSO in a 
given line of sight should have been removed from our re- 
sults. However, there is evidence that the AGN density 
inside clusters is enhanced by a factor of 2 compared to 
non-cluster fields (Molnar et al. 2002; Cappi et al., 2001; 
Sun & Murray 2002), perhaps due to galaxy-galaxy inter- 
actions. Their contribution is not removed by the standard 
background subtraction and thus we need to estimate the 
number of the excess AGN in clusters (compared to blank 
fields). The unobscured AGN are optically identifiable and 
soft X-ray bright. Thus, it is feasible to find them from 
optical catalogs and soft X-ray images, and estimate their 
contribution to HXR, as we describe in Section How- 
ever, the obscured AGN are a difficult problem for the 
HXR studies, because the high obscuration by the torus 
may hide them in the < 10 keV band, and make the op- 
tical detection difficult. At 20 - 80 keV energies, NH has 
no effect and the obscured AGN may give a significant 
contribution in this band. We estimate this contribution 
in Section lOl 

4.1. Unobscured AGN 

A combined ASCA 2-10 keV spectrum of 13 unob- 
scured RIXOS AGNs (Page, 1998) has a photon index of 
1.8±0.1 . Perola et al. (2002) studied nine bright Syl 
galaxies in the 0.1 - 200 keV band using BeppoSAX LEGS 
-I- MEGS -I- PDS data. Their results indicate that in the 
20 - 80 keV band the photon index is 1.8 on average with 
a standard deviation of 0.1, while all best fit values fall 
within the range 1.8±0.2. These two works indicate that 
a slope of 1.8 is a good representation of unobscured AGN 
spectra, and that the extrapolation of the 2-10 keV spec- 
trum up to 80 keV is robust. In our work, when the spec- 
tral information is not available, we use a power-law model 
with a photon index of 1.8±0.2 (at a 90% confidence level) 
as a reference model to estimate the Syl contribution in 
the PDS 20 - 80 keV band data. 

The variable flux level of Syl must be taken into ac- 
count. A study of 113 Syl observed in the ROSAT All 
Sky Survey and in pointed PSPG and HRI observations 
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(Grupe et al. 2001) shows that while a few percent of the 
objects in the sample are transients whose soft band flux 
varies by a factor of 100 in timescales of years, ^ 90% 
of the AGN vary by less than a factor of 2-3. The hard- 
ness ratio analysis is consistent with no spectral variation. 
A study of nine Seyfert 1 light curves in the 2-10 keV 
energy band with RXTE (Markowitz et al. 2001) yields 
results consistent with variability by less than a factor of 
2. Furthermore, they exhibit stronger variability in the 2 
- 4 keV band than in the 7-10 keV band, consistent with 
the ROSAT study. If this trend continues towards higher 
energies, variability by more than a factor of 2 should not 
be common. Thus in our analysis, when simultaneous nor- 
malization level information is not available, we include ± 
50% uncertainty (a factor of 3 variation between lower and 
upper limit at 90% confidence level) for the AGN contri- 
bution to PDS data. 

We searched the SIMBAD database for non-Sy2 AGN 
within 1.3° of the FOV center (we perform a separate 
treatment for Sy2 in Section . We limited the search 
to AGN whose redshifts indicate that they belong to the 
cluster under study. For each cluster, we studied the 3 best 
known objects (see Table 121. We also cross-examined the 
MEGS and PSPG images of the clusters in our sample for 
additional bright point sources. Due to a smaller PSF, we 
used PSPG instead of MEGS to identify point sources, and 
then examined the corresponding MEGS image for excess 
emission in that sky position. We also examined MEGS 
images for additional variable hard band sources, which 
were not visible in PSPG. We assume in the following con- 
servatively that the point sources identified here constitute 
the excess AGN population inside clusters, compared to 
the blank field population, which has not been subtracted 
from the PDS signal. 

The estimation of the AGN contribution to PDS data 
is difficult since there is no spatially resolved hard X-ray 
spectroscopic information for our cluster sample. Thus, 
where possible, we use the MEGS data to obtain the 2 - 
10 keV AGN spectrum and extrapolate it to PDS ener- 
gies. We subtract the local background obtained next to 
the AGN, to ensure similar cluster contributions in both 
source and background data. We include the vignetting ef- 
fect by using ancillary files appropriate for a given off-axis 
angle as provided by the BeppoSAX team. This method 
has the virtue of reducing the uncertainties of the time 
variability of the AGNs. However, due to the wide PSF 
of MEGS, this approach is not accurate for the faintest 
off-axis sources. For those, as well as for the sources out- 
side MEGS FOV, we use PSPG 0.4 - 2.0 keV count rates 
to normalize the reference model, considering the spectral 
and flux variability as described above. The details of the 
AGN modeling in individual cluster fields are given in the 
Appendix. 

4.2. Obscured AGN 

For nearby (z<0.1) Sy2 galaxies to significantly affect 
our results, they need to produce a luminosity of ~ L20-8O 
_ 2943-44 gj.g g-i each PDS pointing (see Section 
I5.1|l . A Ghandra study of the A2104 cluster (Martini et 
al. 2002) revealed five optically unidentified point sources 
(Sy2 galaxies) whose total luminosity reaches 10^'^ erg s~^ 
in the 20 - 80 keV band when using a power-law with pho- 
ton index of 2 to extrapolate from the 2-10 keV band. 



This indicates that Sy2s can in principle affect our results. 

In the unified scheme of AGN and in the X-ray back- 
ground synthesis models it is assumed that the intrinsic 
luminosity distribution of the obscured and unobscured 
objects is the same. Assuming further that the relative 
Syl-to-Sy2 number densities are similar in the field and 
in the cluster environments, we can estimate the 20 - 80 
keV emission of the obscured AGN inside a given cluster 
by multiplying the corresponding Syl contribution esti- 
mated above by a factor of 4. The number of galaxies 
in typical rich clusters is of the order of 100 and thus 
the assumed number of AGN is only a few per cluster, 
which introduces problems of small number statistics. In 
some clusters there are no catalogued Syl and thus no pre- 
dicted Sy2 signal. On the other hand, in clusters A1367, 
A1795, A2029, A2142, A3627, A3667 and Gygnus A the 
Syl based estimate for the non-thermal contribution is 
higher than the observed non-thermal signal. Thus, we 
cannot form a robust Sy2 contamination estimate for each 
cluster, but rather have to resort to a sample average Syl- 
based estimate for the Sy2 contribution. The estimate is 
clearly dominated by Gygnus A and Perseus; excluding 
these sources, the average Sy2 20 - 80 keV band luminos- 
ity of 4 X 10"*^ erg s~^ is similar to the average non-thermal 
luminosity (6 xlO^^ erg s~^) observed in the sample (see 
Fig. Thus, if the assumptions involved in the estima- 
tion are correct, Sy2 galaxies may potentially be a signifi- 
cant source of contamination in the 20-80 keV band. 

The assumption of the similarity of field and cluster 
point source populations can be addressed studying the 
Ghandra analysis of point sources in clusters. These ob- 
servations reach a flux level of 10~^^ erg s~^ cm~^ in the 
2-10 keV band (Martini et al. 2002; Sun & Murray 
2002; Molnar et al. 2002). Following the observation that 
a substantial fraction (~ 50%) of Sy2 are Gompton thick 
(Risaliti et al.), we estimate that the sample average non- 
thermal luminosity yields absorbed 2-10 kev fluxes of 
10^^^ - 10^^^ erg s^^ cm^^, when assuming a power-law 
with aph = 2.0 with NH = 1024-25 ^ni-^ Thus, Ghan- 
dra is sensitive enough to probe a significant fraction of 
the predicted obscured AGN in clusters. However, in sev- 
eral clusters observed with Ghandra (e.g. Molnar et al. 
2002), no such sources were found. Also, the above works 
indicate that the faint point sources in different clusters 
are of different nature, and thus do not support the above 
assumption of substantial field-like Sy2 population in all 
clusters. Gonsequently, the sample average Sy2 contami- 
nation level can only be taken as qualitative. In this work 
we use the quantitative predictions for the fiux of unob- 
scured AGNs, and discuss the possible effects of obscured 
AGNs in the conclusions. 

5. RESULTS 

5.1. Detections 

In order to obtain the count rates of the non-thermal 
emission in 20 - 80 keV band (HXR) we subtracted the 
estimated thermal emission and the unobscured AGN con- 
tribution from the PDS data, and propagated the uncer- 
tainties arising from background fluctuations, PDS data 
statistics, and modeling of the thermal and AGN contri- 
butions (see Fig. Eland TableUI. HXR fluxes vary between 
and 0.1 c/s in 20 - 80 keV PDS band. The statistical 
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uncertainties are similar for different clusters, because the 
PDS signal is dominated by the background and the expo- 
sure times are similar within the sample. The systematic 
uncertainties due to background fluctuations are compa- 
rable to the statistical ones and common for all clusters. 
Thus the uncertainties of the background-subtracted PDS 
count rates are similar in different clusters. The relative 
total uncertainties are quite large, ranging from 10% to 
several 100%. The largest errors correspond to those clus- 
ters with significant Syl contaminations. 

In the sample there are 15 clusters whose 20 - 80 keV 
band signal is not significantly contaminated by Syl (e.g., 
less than 15% of the total signal, thus smaller than the sta- 
tistical errors): A85, A496, A1795, A2029, A2142, A2163, 
A2199, A2256, A3266, A3376, A3562, A3571, Coma, 
Ophiuchus and Virgo. In ~50% of these, the non-thermal 
component is detected at 2 cr level (A2142, A2199, A2256, 
A3376, Coma, Ophiuchus and Virgo). The 4 cr detection of 
the Virgo cluster constitutes a separate case. Virgo is the 
nearest cluster, which renders its data of high S/N, and it 
features the coolest ICM (~2 keV), thus giving the least 
thermal contribution in the PDS band. Furthermore, the 
HXR luminosity of Virgo is one order of magnitude smaller 
than the other 2 a detected clusters, and therefore its hard 
excess is more easily produced by unseen AGNs. We con- 
firm the previously published HXR detections of Coma 
(Fusco-Femiano et al. 1999) and A2256 (Fusco-Femiano 
et al. 2000), albeit at lower confidence level due to the 
level of systematic uncertainties of our work. The HXR 
detection of A2199 would be higher than 2.1 cr, if we as- 
sumed that the steep PSPC spectrum of the AGN in the 
field of A2199 was to be extrapolated to PDS energies, as 
in Kaastra et al. (1999). Again due to our AGN modeling, 
our detection of A3667 is of lesser significance than that 
of Fusco-Femiano et al. (2001). 

All the clusters detected at 2a level, except A2199, ex- 
hibit some degree of merger signatures, i.e. deviations 
from the azimuthally symmetric brightness and tempera- 
ture distributions, reported as follows: A2142 and A3376 
(Markevitch et al. 1998), A2256 (Molendi et al. 2000), 
Coma (Arnaud et al. 2001), Ophiuchus (Watanabe et 
al. 2001) and Virgo (Shibata et al. 2001). The weU es- 
tablished relaxed clusters A1795, A3571 (Markevitch et 
al. 1998), A496 (Markevitch et al., 1999b) and A2029 
(Sarazin et al. 1998) exhibit less significant detections. 
Thus, we divide our sample into two groups: relaxed clus- 
ters (A496, A1795, A2029, A2199, A3571) and merger 
clusters (A85, A1367, A2142, A2163, A2256, A3266, 
A3376, A3562, A3627, A3667, Coma, Ophiuchus), exclud- 
ing Virgo, CygnusA and Perseus, as explained above. As- 
suming that the clusters in both groups lie at the group 
average redshift and that the intrinsic emission models are 
identical inside a group, we formed a weighted mean of 
HXR and its uncertainty for both groups. This yields 0.5 
X 10^2 c g-i and 4.8 X 10-2 c s'^ in PDS 20 - 80 keV 
band for the relaxed and merger group, respectively, i.e. 
the count rate of the merger group is 10 times as high as 
that of the relaxed group. The average redshift of the re- 
laxed group is lower (0.048) than that of the merger group 
(0.058), indicating that the higher count rate of the merger 
group is not due to a distance effect, but rather that there 
are intrinsic differences between the two groups. Due to 
the large systematic uncertainties due to background fluc- 



tuations (1.9-2.7 10^2 c s^^) the detection significance of 
the merger group remains at 2.5ct, while the relaxed group 
count rate is consistent with zero. 

To address the intrinsic emission we assumed that it 
may be modeled with a power-law model with a photon 
index of 2.0, and that only the normalization varies be- 
tween clusters. We normalized this model to the HXR 
values for each cluster to obtain its luminosity I^hxr in 
20 - 80 keV band at the cluster's redshift. Luminosities 
vary from to 10'''' h^Q^ erg s"', most values being in the 
range 10^3-44 ^^-2 gj.g s"' (FigHand TableUJ. Using the 
above model and the average redshift, we converted the 
average count rate of the merger group into luminosity, 
obtaining 8 xlO''^ h^Q^ erg s"', respectively. Thus, the 
data are consistent with a general scenario whereby the 
relaxed clusters have no HXR component, while merger 
clusters do, with a 20 - 80 keV luminosity of ~ io43-44 
Ko erg s-'. 

6. COMBINED SPECTRUM 

Individual cluster signals are of insufficient S/N for the 
purpose of constraining the spectral models. Thus, in or- 
der to obtain information of the average non-thermal clus- 
ter spectrum, we formed an average cluster spectrum by 
co-adding bin-by-bin the PDS counts of each cluster whose 
20 - 80 keV band signal is contaminated by less than 10% 
by Syl (see EH), i.e. A85, A496, A1795, A2029, A2142, 
A2163, A2199, A2256, A3266, A3376, A3562, A3571, 
Coma and Ophiuchus. The co-added exposure time is 560 
ks. To avoid artificial overestimation of the uncertainties, 
we did not propagate the uncertainties of the individual 
spectra, but rather used the combined spectrum to deter- 
mine the Poissonian uncertainties. The level of systematic 
uncertainty due to background fluctuations in the 20-80 
keV band is ^15% of the background-subtracted PDS sig- 
nal. Since we have no information on the energy depen- 
dence of this quantity, we assumed it to be a constant 15% 
in the 12 - 115 keV band. Combining this with the uncer- 
tainty of AGN contamination, we arrive at 20% systematic 
uncertainty, which we use in the following analysis. 

In an attempt to account for the total thermal contri- 
bution, we first fitted the 12 - 20 keV band data with a 
mekal model keeping metal abundance at 0.3 Solar and the 
redshift at the median of 0.06. In this band, the typical 
NH of 10^° cm-2 has no effect, and thus we exclude the 
absorption from the model. The best fit temperature is 
consistent with the median of 7.8 keV of the sample, im- 
plying that the non-thermal emission does not dominate 
in the 12 - 20 keV band. The non-thermal excess on top 
of the thermal model is clearly evident: at 100 keV, the 
thermal model underpredicts the signal by 4 orders of mag- 
nitude (see FiglSlfor the thermal contribution in the final 
best fit model). Fitting the 12 - 115 keV band data with 
only a mekal model, we obtain a statistically unacceptable 
fit with unrealistically high temperature of 26 keV. This 
further confirms the existence of an additional hard X-ray 
component. 

We introduced a power-law component to the 12 - 115 
keV band fit, allowing the photon index and the normal- 
ization to vary, together with the mekal temperature and 
normalization. The best-fit is formally acceptable, with 
xVdof of 10.9/11, yielding a photon index of 2.8toi (Fig. 
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7). The typical AGN photon index of 1.8 is ruled out at 
98% confidence level, which argues against significant Sy2 
contamination in the hard X-ray signal. 

However, the non- thermal emission in this model in the 
12 - 20 keV band is high, ^ 50% of the total. In order 
to stiidy the relative contribution of the thermal and non- 
thermal components, we examined the central 8' MEGS 
data of the largest contributors to the thermal emission 
in the sample, i.e. Goma and Ophiuchus by fitting the 
MEGS 2-10 keV data with a mekal -I- a power-law with 
aph fixed to 2.8. The allowed Icr upper level for the non- 
thermal contribution, extrapolated to 12 - 20 keV band, 
is below 1%. Thus, the above best fit model requires that 
the non-thermal emission is extended and negligible in the 
central 8'. To confirm that this is the case, one needs to 
perform spatially resolved hard X-ray spectroscopic anal- 
ysis on the cluster sample, which is currently not possible. 

Alternately, if we assume that most of the HXR origi- 
nates from cluster centers, the MEGS data require a harder 
spectrum for it: Fixing aph to siiialk^r values, and keep- 
ing NH at Galactic values, the mekal -|- power-law fit to 
2-10 keV MEGS data of Goma and Ophiuchus allow big- 
ger contribution from the non-thermal model in the 12 - 
20 keV band. Also, forcing the non-thermal component in 
the PDS 12 - 115 keV fit to be harder, decreases the non- 
thermal contribution in 12-20 keV PDS band and with 
ctph < 1.5 the non-thermal flux at 12-20 keV in the best- 
fit PDS model is below 10% of the thermal, consistent with 
the MEGS data of Goma and Ophiuchus. On the other 
hand, the decreasing non-thermal contribution in the 12 
20 keV band requires higher temperatures for the best-fit 
PDS model, and at aph < 1.3 it exceeds the highest tem- 
perature of the cluster sample. Thus, assuming that the 
non-thermal emission comes mainly from the cluster cen- 
ters, its photon index is limited within 1.3-1.5. However, 
such hard slopes yield poor fits to the PDS data. Keeping 
aph = 1.3 the model has x^/dof = 20.9/12 and systemat- 
ically excecids the; data by 20-40% above 70 keV energies. 
Note that assuming an obscured AGN model with aph 
= 1.8 and NH = 10^^ cm~^ would yield consistence be- 
tween MEGS and PDS, since the high absorption in 2-10 
keV band would hide the non-thermal contribution from 
MEGS data, but would allow it to dominate in PDS band. 
However, this forced fit is also significantly worse (x^/dof 
= 17.2/12) than that of the free aph fit. Thus the data 
indicate that the steeper slope (2.4 3.1) and thus the ex- 
tended distribution of the non-thermal emission are more 
likely. The Sy2 galaxies inside clusters are concentrated 
in the central high galaxy density regions and thus the in- 
dication for extended nature of the non-thermal emission 
also argues against the Sy2 origin of the PDS signal. 

7. MODELS AND DISCUSSION 

Most models for the HXR emission require acceleration 
of cluster electrons to supra-thermal and/or relativistic ve- 
locities. Large-scale acceleration is natiirally provided by 
merger shocks, and our findings (higher HXR detection 
significance of the merger clusters compared to relaxed 
clusters, see Section 5) arc consistent with this basic as- 
sumption. A strong merger accelerates electrons to rela- 
tivistic velocities and consequently the Gosmic Microwave 
Background photons may experience Inverse Gompton 
scattering (IG/GMB) from these electrons, thus producing 



hard X-ray emission. Within the framework of the merger 
acceleration, the observed photon index (2.4 3.1) of the 
combined PDS spectrum in this work implies a power-law 
form for the differential momentum spectra of the relativis- 
tic (^ GeV) electrons with a slope of /i = 3.8 - 5.2. Right 
after the first acceleration event the primary electron dis- 
tribution is predicted to be harder (^ 2-2.5, see Miniati et 
al., 2001) but the electrons loose energy rapidly and their 
spectrum in the GeV range steepens into consistence with 
that derived from the observed PDS spectrum. 

The IG/GMB model requires a confinement of the rel- 
ativistic electrons in clusters, which can be achieved by 
the cluster magnetic fields. In the presence of magnetic 
fields, the relativistic electrons produce synchrotron emis- 
sion at radio wavelengths. Thus, the model naturally pre- 
dicts a connection between the non-thermal hard X-rays 
and radio emission. We have indicated above a connec- 
tion between the non-thermal hard X-ray emission and 
cluster mergers, which in turn predicts a connection be- 
tween cluster mergers and radio emission. Indeed, clusters 
with a large scale (> 1 Mpc) radio halo possess merger sig- 
natures such as substructure in the X-ray brightness and 
temperature distribution and absence of cooling flows (see 
Feretti 2003 and references therein). Also, the diffuse radio 
emission is more common in clusters with higher X-ray lu- 
minosities (Giovannini et al. 1999), perhaps due to energy 
input by recent mergers, as in hydrodynamic simulations 
(Sarazin et al. 2002). Thus, the observed connections of 
cluster mergers with radio emission and with non-thermal 
hard X-rays support the IG/GMB scenario whereby merg- 
ers provide higher temperatures and luminosities as well as 
stronger shock acceleration, and thus stronger radio and 
non-thermal X-ray emission. Within this framework, the 
observed spectral index of the combined PDS spectrum in 
this work (1.4 - 2.1) equals that of the radio spectral in- 
dex of the synchrotron spectra. Indeed, most radio-halo 
cluster observations typically yield spectra with indexes 
in this range (e.g. Feretti et al. 2001; Fusco-Femiano et 
al. 1999), further strengthening the case for IG/GMB. In a 
forthcoming paper we will examine the connection between 
the HXR and the radio emission in a sample of clusters in 
more detail. 

The case for IG/GMB is challenged by the clusters 
A2163, A3266 and A3562 with merger signatures featuring 
less significant HXR detections from the rest of the merger 
group. Also, HXR for a relaxed cluster A2199 is detected 
with 2a confidence. This implies that the merger is not the 
only factor responsible for the non-thermal emission in all 
clusters, which gives room for other models. Even though 
there is a possibility of Sy2 contribution in the PDS signal 
(see above), the co-added spectrum is steeper than those 
observed in AGN and the indicated extended distribution 
of HXR is also contrary to the central concentration of 
AGN in clusters. Thus the current data argue against sig- 
niflcant contamination by obscured AGN in our sample. 
The non-thermal bremsstrahlung model (e.g.. Sarazin & 
Kempner 2000, Dogiel 2000) predicts spectral slopes of 
HXR consistent with our observations, and thus cannot 
be ruled out by the flt to the current data. However, 
bremsstrahlung is a very inefficient process (see, e.g., Pet- 
rosian 2001, Timokhin et al. 2003), and the huge amount 
of energy input needed to produce the observed level of 
hard X-ray emission is ruled out in cases like Goma by 
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X-ray observations (Petrosian 2001). 

The secondary electron models usually predict a harder 
[oiph ~ 1.5-1.75) IC/CMB spectrum (Colafrancesco & 
Blasi 1998; Blasi & Colafrancesco 1999; Miniati et al. 
2001) than that indicated by the present observations. If 
we assume that most of the non-thermal hard X-ray emis- 
sion originates from the central regions of the clusters, its 
spectrum is required to be hard {aph = 1.3-1.5, see Sec- 
tional). In this case, the spatially concentrated hard X-ray 
emission is consistent with the secondary models, which in- 
volve the production of secondary electrons via collisions of 
relativistic protons which are bound to the cluster gravita- 
tional potential wells (Colafrancesco & Blasi 1998). How- 
ever, this model does not fit well the PDS data at highest 
energies. 

8. CONCLUSIONS 

We have studied the hard X-ray emission in 20 - 80 
keV energy band in a sample of clusters using the Bcp- 
poSAX PDS instrument. After removing the contributions 
from the cluster thermal component and from unobscured 
AGN, in ~50% of the mildly AGN-contaminated clusters 
the non-thermal component is detected at 2 cr level, the 
clusters being A2142, A2199, A2256, A3376, Coma, Ophi- 
uchus and Virgo. All the clusters detected at 2a level 
exhibit some degree of merger signatures, i.e. deviations 
from the azimuthally symmetric brightness and tempera- 
ture distributions (except for A2199). Averaging the PDS 
20 - 80 keV count rates of the relaxed and merger clusters 
obtains a 2.5cr detection for the merger group, while the 
relaxed group count rate is consistent with zero. Assum- 
ing a power-law emission model with a photon index of 
2.0 at the group average redshifts, the average count rates 
are consistent with a scenario whereby the relaxed clusters 
have no HXR component, while mergers do, with a 20 - 
80 keV luminosity of ~ 10*3-44 j^-a gj.g g-i^ 

The co-added spectrum of our sample yields a best-fit 
photon index of 2.8^0 4 for the non-thermal emission in 
12 - 115 keV band, and we find indication that it has ex- 
tended distribution. These indications argue against sig- 
nificant contamination from obscured AGN, which have 
harder spectra and centrally concentrated distribution. 



The indicated connection between cluster mergers and 
the non-thermal hard X-ray emission is consistent with 
the Inverse Compton scattering of the Cosmic Microwave 
Background photons with merger-accelerated population 
of relativistic electrons. In this framework, the observed 
photon index is consistent with a scenario in which a strong 
acceleration event and consequent strong IC/CMB energy 
losses take place. In this scenario the measured hard X- 
ray slope corresponds to a differential momentum spectra 
of the relativistic electrons with a slope of yit = 3.8 - 5.0. 
The consequent synchrotron emission spectrum expected 
from the same electron population has a spectral index of 
1.4 - 2.1, consistent with radio halo observations of many 
merger clusters. 

The observed slope of the HXR spectrum is also 
consistent with the predictions of the non-thermal 
bremsstrahlung model. Even though this fit cannot be 
ruled out by the current data, the bremsstrahlung model 
seems to face a strong energetics problem which does not 
make it a viable physical scenario. 

Assuming that most of the non-thermal signal originates 
in the central regions of clusters, the HXR spectrum is 
forced to be harder, with a slope 1.3-1.5, which turns 
out to be consistent with secondary electron models. How- 
ever, this model provides a worse fit to PDS data and is 
thus disfavored by the statistical fit over the primary elec- 
tron IC/CMB model. 

In conclusion, spatially resolved hard X-ray spec- 
troscopy is needed to disentangle between primary and sec- 
ondary electron models for non-thermal hard X-ray emis- 
sion in clusters of galaxies. 
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Fig. 1. The PDS count rate distribution of the selected 164 fields in 20 80 keV band. The histogram shows the data 
when using the standard method of subtracting the background of both off-sets. The solid curve is a sum of gaussians 
which represent the data points with the value of the data point as the centroid and the uncertainty as a. The dash-dot 
and dashed lines show the corresponding sum of gaussians, when using only positive or negative off-sets for background 
subtraction. 
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Fig. 2. — The 20-80 keV count rate difference between positive and negative offset background pointings for eacli cluster 
exposure. Tlie daslied line shows the average difference. 
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Fig. 3. — The non-thermal signal and la uncertainties in PDS 20 - 80 keV band after subtraction of the contributions 
from the background, thermal gas and AGN in the field, and after propagating uncertainties due to these subtractions. 
The dotted vertical line separates the relaxed clusters (left) from the rest (right). 
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Fig. 4. — The luminositites of the non-thermal emission in PDS 20 - 80 keV band at la confidence level, obtained by 
using a power-law model with a photon index of 2.0. The dotted and dashed lines show the allowed 90% upper limit for 
HXR luminosity in the relaxed and merger clusters. Note that A3571 and A2029 values are negative at upper la level 
and thus excluded from the logarithmic plot. 
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Fig. 5. — The combined spectrum of all the clusters not significantly affected by AGN. The lines show the unfolded model 
components while the crosses show the data and 1 a errors (including 20% systematics). The solid line shows the total 
model. The dotted line shows the thermal contribution. The dashed line shows the best fit power-law of a^h — 2.8. 
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Table 1 
Results for HXR 



name 


PDS CR" 


a 


thermal" 


AGN"* 


HXR^ 


CLhxr 
a 


HXRqo^ 










2 7 


4.1±0.8 


- 


7 6+"^-^ 

' •'J-4.4 


1 7 

J. . 1 


14.8 


If) 7+6-3 
J^'-'- ' -6.3 


21.1 




-L .LI— 1— O. 1 














A496 


3.7±3.4 


1.1 


1.5±0.4 




2.3±3.5 


0.7 


8.0 


1.3±2.0 


4.6 


A1367 


2.1±3.1 


0.7 


0.4±0.2 


1 0+1^1 


0.7±3.1 


0.2 


5.9 


0.2±0.8 


1.6 


A1750 


1.6±4.0 (-) 


0.4 














A1795 


3.0±2.9 


1.0 


3 7+nI 

' -0.6 


0.3±0.1 


-1.0±3.0 




3.9 


-2.0±6.0 


7.8 


A2029 


2.0±4.8 


0.4 


8.9±1.3 


0.2±0.1 


-7.1±5.0 




1.0 


-22.0±15.3 


3.2 


A2142 


25.0±4.4 (+) 


5.7 


11.9+J-l 

— 1.8 


3.0+^f 

"^—0.6 


10.1±4.8 


2.1 


18.0 


42.0±19.7 


74.6 


A2163 


8.3±3.3 


2.6 


6.5=b0.9 




1.8±3.4 


0.5 


7.4 


38.6±73.1 


159.4 


A2199 


9.2±3.4 


2.7 


1.6=b0.4 


'-'•"-0.3 


7.1±3.4 


2.1 


12.7 


3.4±1.6 


6.0 






4 




7 c_|_9 q 


9 


1 9 


J-O.unztj.o 


99 






2.1 


















y.ozto.o 




0.4zti.U 


U-^-O.l 


z.yito.y 


U.o 


y.o 


4.oitD.4 


10. 4 


A3376 


9.5±3.0 


3.2 


0.4±0.2 


0.9+°-3 

— u. z 


8.2±3.0 


2.7 


13.1 


9.2±3.4 


14.7 


A3562 


6.8±5.1 (-) 


1.3 


0.9±0.6 


5.8±5.2 


1.1 


14.3 


8.0±7.1 


19.7 


A3571 


-1.7±4.6 (+) 




S.O'^ 




-4.8±4.6 




2.8 


-4.1±3.9 


2.4 


A3627 


11.4±3.4 


3.3 


6.8±0.6 




1 c+3.8 
i-0_4.2 




4.7 


2+°-*^ 

"•^-0.7 


0.8 


A3667 


8.2±2.6 


3.2 


4.4±0.4 


1 7+0.6 
' -0.4 


2.1±2.7 


0.8 


6.4 


3.2±4.1 


9.9 


Coma 


40.2±3.4 


11.8 


30.9±0.8 


0.9±0.6 


8.4±3.5 


2.4 


14.3 


2.6±1.1 


4.4 


CygA 


58.6±3.4 


17.1 
















Ophiu 


75.0±3.9 


19.3 


66.2^2° 






2.0 


15.5 


3.7±1.8 


6.6 


Perseus 


54.7±3.5 


15.7 


37.6±d.6 












PKS07 


3.5±3.3 


1.1 
















RXJOl 


0.8±3.0 


0.3 
















RXJ13 


6.8±4.4 (-) 


1.5 
















Virgo 


27.5±5.5 


5.0 


0.3±0.2 


4.5±0.7 


22.8±5.6 


4.1 


31.9 


0.17±0.04 


0.24 


z3146 


1.4±4.4 


0.3 

















^^10-2 c s-1 in PDS 20 - 80 keV band 
^lO'*^ hgtf erg s-i in 20 - 80 keV band 
'^predicted 



Note. — The count rates (PDS CR) are obtained with PDS in 20 - 80 keV band using standard background subtrac- 
tion method (default) or using only pos (+) or ncg (-) background pointing and correcting for the systematic effect. 
The errors include both statistical and systematic uncertainties at 1 a level. CL^et gives the confidence level of source 
detection, "thermal" gives the thermal model prediction in 20 - 80 keV band, normalized to PDS 12 - 20 keV data, 
together with 1 a uncertainty due to PDS photon statistics in the 20 - 80 keV band. "AGN" gives the estimated AGN 
contribution to PDS 20 - 80 keV band. HXR and CLhxr give the non-thermal, AGN subtracted count rate with 1 a 
errors, and its confidence level. HXRgo gives the 90% confidence upper limit of the count rate of the non-thermal emis- 
sion. Lhxr gives the luminosity of the AGN subtracted non-thermal component, obtained by normalizing a power-law 
with photon index aph = 2.0 to HXR, and its la uncertainties. Lhxr90 gives the the 90% confidence upper limit of 
the HXR luminosity. 



Table 2 
Thermal models 



name 




z 


Nh 


I1/I2 


['] 


rc2 
['] 






T 

[keV] 


ab 
[solar] 


Refs. 

/3,T,L 


A85 


0, 


.052 


2.7 


0.08 


3.9 


0.72 


0, 


.66 


10.3 


^ /-.-i-r) 4 






1,8 ,20 


A348 


0, 


.274 


3.1 




STANDARD 






3.7 








7,12,12 


A496 


0, 


.033 


4.2 


0.05 


4.0 


0.55 


0, 


.65 


4.3 


*-'-0.2 






1,10,20 


A1367 


0. 


.021 


2.4 




10.0 




0, 


.61 


3.7 


r» ^,-1-1-0 10 

3.69^°;}° 






1,9 ,20 


A1750 


0, 


.086 


2.4 




STANDARD 






3.2 


A A jr,-\-f\ OA 

4.46lg-Jl 


0.30 


+0.10 
-0.10 


7,9 ,13 


A1795 


0, 


.062 


1.0* 


0.07 


3.5 


0.82 


0, 


.79 


12.1 


7.8+};° 






1,8 ,20 


A2029 


0, 


.077 


3.2* 


0.05 


2.8 


0.68 


0, 


.71 


17.2 


-1-1 n 

9-iti;o 






1,8 ,20 


A2142 


0, 


.089 


4.1* 


0.07 


4.8 


1.2 


0, 


.79 


22.5 








1,8 ,20 


Azibo 


ri 

V. 


.zU3 


1 1 n 




1.6 




0, 


.73 


/too 
43.3 


11 c 
ii.O 






oil OA 

2,ii,2U 


A o 1 f\r\ 
Azi99 


U, 


.030 


U.OD 


A 1 

O.lo 


3.2 


0.81 


0, 


.66 


A A 

4.4 








1 1 A OA 

1,10,20 


A2256 


0, 


.058 


4.5* 




5.3 




0, 


.83 


8.2 


a n'7+0.12 
6.97I0.12 


0.26 


+0.02 
-0.02 


1,9 ,20 


A2390 


0, 


.228 


7.0 




0.47 




0, 


.60 


26.3 


A 0+0.8 

9-8-0.7 


n 3+°-^ 


3,14,20 


A32dd 


0, 


.055 


1.6* 




5.7 




0, 


.74 


8.1 


„ „„+0.30 

8-97Io.3o 


0.22 


+0.03 
-0.03 


1,9 ,20 


A3376 


0, 


.046 


4.4* 




STANDARD 






3.0 


AA + 0.13 
-\ MM 


0.23 


+0.04 
-0.04 


7,9 ,20 


A3562 


0, 


.050 


3.8 




1.2 




0, 


.47 


7.1 


r 1 +0.3 

^•J^-0.3 


0.39 


+0.08 
-0.08 


1,15,20 


A3571 


0, 


.040 


4.1* 




2.6 




0, 


.61 


9.7 


6.9+^'2 






1,8 ,20 


A3627 


0, 


.016 


21.9 




10.0 




0, 


.56 


3.8 


6.28t°;l« 


0.27 


+0.02 
-0.02 


4,9 ,20 


A3667 


0, 


.053 


4.8 




3.1 




0, 


.54 


14.4 








1,8 ,20 


Coma 


0. 


.023 


0.9 




10.1 




0, 


.71 


8.8 


01 +0-16 


0.21 


0.03 
-0.03 


1,18,20 


Cygnus A 


0, 


.057 


36.1 




0.17 




0, 


.47 


15.5 


6.9t}j 


0.67 


+0.12 
-0.10 


1,8 ,8 


Ophiuchus 


0, 


.028 


20.3* 


0.65 


5.8 


1.7 


0, 


.71 


13.5 


9.itS;? 


0.49 


+0.08 
-0.08 


1,13,13 


Perseus 


0, 


,018 


14.8 


0.02 


13.1 


2.0 


0, 


.75 


9.6'' 


0.00_o 18 


0.41 


+0.02 
-0.02 


1,19,13 


PKS0745-191 


0, 


.103 


42.4 




STANDARD 








8 5+° '' 


0.38 


+0.03 
-0.03 


7,17,- 


RXJ0152.7-135.7 


0, 


,831 


1.6 




STANDARD 






11.6 


fi r,+2.9 
D.O_2.0 


0.5l 


0..5 
0.4 


7,16,16 


RXJ1347.5-1145 


0, 


.451 


4.8 




0.14 




0, 


.56 


58.4 


14.3+1-^ 


0.5l 


0.2 
0.2 


5,14,14 


Virgo 


0, 


.0036 


2.5 




2.2 




0, 


.45 


0.7'^ 


9 qr+0:b6 
z.oo_o.og 


0.49 


+0.06 
-0.06 


6,13,13 


z3146 


0. 


.291 


3.0 




STANDARD 






28.3 


7 q+0.9 

'•"5-0.8 


n 3+°-^ 


7,14,14 



'^lO^^ hgjf erg s-i in 0.1 - 2.4 keV band 
''power-law component removed 

Note. — The Nnvalues [ 10^" atoms cm"^ ] are based on Dickey & Lockman, except for the ones marked 
with * which are taken from fine beam HII survey (thin filter) of Murphy et al. (in prep.) The /3 model 

parameter reference (7) corresponds to STANDARD model of /? = 2/3 and rd = 0.2 Yi^^ Mpc, due to lack 
of proper reference. The abundances, where marked, are taken from the temperature references, otherwise 
0.3 Solar is assumed. The unabsorbed hmiinosities Lyg' [lO^** h^g^ erg s"'^] in 0.1 - 2.4 keV band are obtained 
by using the (5 models to extrapolate the luminosities taken from the references papers. 

References.— (1) Mohr et al. 1999; (2) Vikhlinin et al. 1999; (3) Bohringer ct al., 1998; (4) Bohringer 
et al. 1996; (5) Schindler ct al., 1997; (6) Bohringer et al. 1994; (7) STANDARD; (8) Markcvitch ct al. 
1998; (9) dcGrandict al. 2001; (10) Markcvitch ct al. 1999; (11) Markcvitch ct al. 1996; (12) Colafrancesco 
et al. 2001; (13) this work, (14) Ettori ct al. 2001; (15) Ettori et al. 2000; (16) Delia Ceca et al. 2000; 
(17) de Grandi et al. 1999; (18) Hughes et al. 1993; (19) Allen et al. 1992; (20) Ebeling et al. 1996 
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Table 3 

Observation log and AGN information 
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APPENDIX 

THE DETAILS OF THE MODELING OF THE UNOBSCURED AGN IN INDIVIDUAL CLUSTERS 

A85, A496, A2163,A2256: There are no catalogued Syl within these clusters, nor strong point sources in MECS or 
PSPC in the field. 

A1367: MECS spectrum of AGN NGC 3862 provides constraints on the reference model, yielding a PDS estimate of 
1.0^03 10~^ c s~^. This value is consistent with HXR and thus the HXR estimate will be very uncertain. The PSPC count 
rate of NGC 3884 is negligible (< 1 %) of that of NGC 3862 and thus will not change the AGN contribution estimate. 

A1795: There is a Seyfert I galaxy IE1346+26.7 5' off-axis. Normalizing the reference model to MECS data gives 
a PDS estimate of 0.3^0,1 10"^ c s"^ In the cluster center there is LINER PKS 1346+26. Its flux estimate cannot be 
given due to the projection with the bright cluster center. However, LINERS usually have 2-10 keV luminosities 1-3 
orders of magnitude smaller than classical Seyferts (Terashima et al., 2002). AGN RXJ1343.9+2712 is outside the PSPC 
image and thus an flux estimate cannot be given. The fleld of AI795, together with Coma, is unusual in its large number 
of AGN/QSO. For AI795 this is probably not a problem though, because HXR is negative (consistent with 0) and likely 
not significantly contaminated by any AGN. We will use the IEI346+26.7 estimate in the following. 

A2029: The MECS data of an AGN QSO J1511+057 at 8' off-axis constrain the reference model yielding PDS estimate 
of 0.2to\ 10-2 c s-^ A Syl JVAS B1509-F054 is located at 29' off-axis and thus outside MECS FOV. It is undetected 
in the PSPC, and the upper limit of statistical uncertainties allow 10% of the PDS estimate of QSO J1511-I-057 which is 
negligible. 

A2142: A Seyfert 1 galaxy IC 1166 is outside the PSPC field and thus the flux estimate cannot be given. 4 ' ofi^-axis 

from the cluster centre there is a Seyfert 1 galaxy IE 1556+27.4 whose photon index is apt = 1.9 as observed with ASCA 
(Markevitch et al., 1998). The MECS data gives consistently aph = 1.8±0.1. Using this and including the statistical 
uncertainties of the MECS data, the extrapolated PDS 20-80 keV count rate is 2.2tol 10"^ c s'^, consistent with PSPC 
estimate of 2.9 10"^ c/s. In A2142 there is another Seyfert I galaxy, QSO BI557+272. At 17' off-axis, the source is 
quite diffuse in MECS, but the data extracted from a 4' circle around the source still gives adequate constraints on the 
normalization of the power-law model. Including statistical imccrtaintics, the PDS prediction is O.S^q ^ 10^^ c s~^. Thus 
the combined AGN contribution to PDS 20 - 80 keV band is 3.0lg;^ lO'^ c s^\ or 30% of the HXR signal. 

A2199: There are no catalogued Syl in A2199, but X-ray imaging reveals a bright Seyfert I galaxy, RXS J16290+4007 
of rcdshift 0.3, located 35' off-axis and thus outside the MECS FOV. Its PSPC spectrum is not consistent with the 
reference model. The data can be modeled with a combination of thermal and a power-law model, yielding a temperature 
of 0.1 keV and a photon index of 2.4+0.1, which gives negligible contribution to HXR. Allowing for spectral hardening 
towards higher energies, we fitted the PSPC data with our reference model + mekal. The fit is bad. but however yields 
an estimate for HXR of 0.61q'3 10^^ c s~^, 8% of the HXR, thereby decreasing the detection confidence slightly. To 
be consistent with the treatment of Seyfert 1 galaxies in the rest of the sample, we assume the harder spectrum in the 
following. 

A3266: A Syl J043829.3-614759 at the edge of PSPC gives constraint to the reference model, yielding a PDS estimate 

of 0.2tol 10"^ c s-i, 5% of the HXR. AGN C3266-I2 and E3266-3 are outside MECS FOV and undetected in the PSPC. 
The upper limit allowed by statistical imccrtaintics of PSPC data, fitted with the reference model, yields a PDS estimate 
below 1% level of HXR and thus negligible. In the PSPC there is a bright point source IRXS J043356. 7-6 12909 at 
04'*33""56.70*, -61°29' 09.5". However it is not visible in MECS, implying that the source is either very soft or variable, 
and very faint during BcppoSAX observation. Either way it gives no contribution to PDS. Thus we keep the estimate of 
J043829.3-6I4759. 

A3376 In the field there is a bright point source IRXJ J060II3. 0-401643 at 06''0I'"32", -40°16'55.7", 18' off-axis. There 
is no information available on its nature. MECS constraints on the spectral slope are poor, but the data however give 
good constrain on the normalization of the power-law model, when fixing the photon index to 1.8+0.2. With this model, 
we obtain PDS estimate of 0.3+0.1 10"^ c s"^. A bright source at off-axis 24' co-incides with QSO IWGA J0600.5-3937 
and source PKS 0558-396. Using MECS spectrum we obtain PDS estimate of O.Gtol 10"^ c s^^. A Syl J055850.3-403848 
within A3376 at off-axis of 50' has count rate much below the above sources and thus has no effect on the combined 
estimate of 0.9lg-^ 10"^ c s-\ ~10 % of the HXR. 

A3562: Syl I'EI335.I-3I28 at 60' off-axis is undetected in PSPC, implying a negligible HXR contribution. A poor 
cluster SCI329-3I3 is included in the PDS FOV. 

A3571: There are no catalogued AGN or QSO in the cluster. There is a bright point source HD 119756, an X-ray 
binary HD 119756 in the field at RA, DEC = 13'*45'"41.5% -33°02'32". It is obscured by the MECS calibration source 
and PSPC spectrum indicates thermal spectrum with T = 0.4 keV with no evidence for power-law component. Thus for 
A357I we estimate negligible PDS contribution from point sources. 

A3627: Close to the edge of MECS there is a projected Seyfert I galaxy IWGA J16II.8-6037. Using the MECS 
spectrum we normalized the reference model and obtained PDS estimate of 6.ll^ g 10"^ c s~^. This is consistent with 
the HXR estimate which will thus be very uncertain. 

A3667: Seyfert 1 galaxy FRL 339 within A3667 is close to the edge of MECS FOV. The MECS spectrum provides 
constraint for the power-law component as I.9lg'2, consistent with the reference model. Using the MECS data we 
normalized the reference model and obtained PDS estimate of I.OIq 2 10"^ c s~^. In the PSPC image there are two other 
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bright non-catalogued point sources, 2E 2007.4-5653 at RA, DEC = 20'' 11" 28.6", -56°44'13" and IRXS J201455.6- 
565833 at RA, DEC = 20'n4™55.6% -56°58'33". The former is undetected in MECS because it is projected at the bright 
cluster center and the latter is outside MECS FOV. Both are classified as X-ray sources. The PSPC data of 2E 2007.4- 
5653 is consistent with the reference model, and gives a PDS estimate of 0.7+^J 10"^ c s-\ The PSPC data of IRXS 
J201455. 6-565833 is not consistent with the reference model, and requires a steeper photon index. The PDS prediction 
with this model is insignificant. The combined AGN contribution (using FRL 339 and 2E 2007.4-5653) is 1.7^^ j 10"^ ^ 
s~^, ~ 80% of HXR. The point source contamination is not discussed in the report on the marginal hard excess of A3667 
(Fusco-Femiano et al. 2001). 

Coma: The well known Seyfert 1 galaxy X-Comae is just at the edge of the MECS FOV. Fusco-Femiano et al. (1999) 
used MECS data to show that the allowed upper flux level for X-Comae is ^ 15 % of the hard X-ray excess component at 
2-10 keV, using a power-law component with Uph = 1-8. We used PSPC data to check the normalization of the reference 
model. Including the spectral and flux level variation uncertainties we obtain a PDS 20 - 80 keV estimate of 0.9to'6 
10-2 c/s or 10 % of the 20-80 keV HXR emission, consistent with Fusco-Femiano et al. (1999). AGN IE 1258+28.9 at 
off-axis of 50' is obscured by the PSPC mirror support structure. The useful data still indicates a similar count rate as 
for X-Comae, indicating a significant contribution to PDS. However, the nature of the source is not well known and the 
extrapolation towards higher energies is not justified. Syl J125710. 6+272418 is imdetected and thus its contribution is 
negligible compared to that of X-Comae. Coma field contains an unusually high number, 26, of catalogued AGN/QSO, 
perhaps because Coma field is better studied than others. However, to make up all the HXR, seven objects like X-Comae 
are needed and this is ruled out for Coma based on the PSPC image. 

Cygnus A: There is a powerful radio galaxy QSO B1957+405 in the center of Cygnus A. We extracted central 2' MECS 
spectrum and modeled it as a sum of mekal and self-absorbed power-law, both absorbed by the galactic Njj. The best fit 
photon index aph — 1.9lg'2 consistent with the ASCA result (Markevitch et al. 1998), and with our reference model. 
We thus used the 2' MECS data to normalize the reference power-law model, including a mekal model with temperature, 
metal abundance and normalization as free parameters. Extrapolating the resulting power-law model to higher energies 
we obtained the PDS estimate of IQ-^ c s"^. This is consistent with the total observed PDS emission in this 

band and thus the HXR estimate will be uselessly uncertain. We thus reject Cygnus A from further analysis. 

Ophiuchus: There are no catalogued AGN or QSO in Ophiuchus. In PSPC image there is a bright point source RXS 
J171209.5-231005 at RA,DEC = 17''12™09^ -23°09'50", classified as X-ray source. It is undetected in MECS due to 
projected bright cluster center in the line of sight. The PSPC spectrum exhibits a 2 component spectrum, consisting of 
a thermal one with T ~ 1 keV and a very steep (a > 3) power-law, which gives negligible contribution to PDS HXR. 

Perseus: Perseus hosts a well known AGN NGC1275 in the center. HEAO I observations revealed a non-thermal 
component in Perseus data at 20 - 50 keV band (Primini et al. 1981). The excess was modeled with a power-law model 
whose best fit photon index Uph is 1.9±0.3 at 90% confidence. They also report that the source exhibits no significant 
variations above 25 keV in time scale of 4 years. PDS data are of high enough quality to perform two component fit, if 
we fix photon index (= 1.9) and mekal abundance (= 0.3). The resulting temperature 6.3±0.4 keV is identical with the 
Ginga value (Allen et al. 1992). The power-law component has 25 - 40 keV luminosity of 1.8±0.5 10*^ erg s~^, 4 times 
smaller than The HEAO I value, implying variability on a time scale of 20 years. 

Due to the high brightness of cluster thermal emission in the center, compared to that of NGC1275, the central 2' 
MECS data do not provide decent constraints on the internal Nh or the power-law slope of the AGN. However, modeling 
the central MECS data with mekal + the above power-law component reveals that the power-law model given by PDS 
data contributes only a few % of the total emission. This component modifies the total model only slightly and the fit is 
acceptable. If let free, the allowed upper limit for the normalization of the power-law component is 3 times as high as the 
best value given by PDS data. Thus, the data are consistent with all of the non-thermal emission coming from NGC1275. 
We thus reject Perseus from further analysis. 

Virgo: Virgo has an active nucleus M87 and a jet in the center. XMM-Newton data yields power-law slopes of 2.2±0.2 
and 2.5±0.4 for the nucleus and the bright knot in the jet (Bohringer et al. 2001) at 90% confidence, and no indication 
of excess absorption. The 2' MECS data do not provide good constraint on the slope of the non-thermal component. We 
thus fit the 2' MECS data with a model consisting of mekal and power-law, both absorbed by the galatic Nh , fixing the 
photon index to 2.3, based on the XMM-Newton observations (Bohringer et al. 2001), thus obtaining the normalization 
and its uncertainty for the power-law component. We determined the thermal component of Virgo by using the above 
determined central power-law model together with a mekal when fitting 0-8' keV MECS data. The best fit parameters T = 
2.35+0.04 keV and abundance 0.49+0.04 Solar are consistent with XMM-Newton results. Letting only the thermal model 
normalization as a free parameter, we then normalized this model to PDS FOV using 12-20 keV PDS data. According 
to this model, M87 contributes 17+3% of the non-thermal emission in 20 - 80 keV band. Allowing spectral variability 
for M87, we repeated the above exercise keeping Uph at 2.0 and 1.7. The resulting M87 contribution to HXR is 20 - 30% 
and 30 - 50%, respectively. Unless the spectrum of M87 has a strong hard excess, the non-thermal PDS signal of Virgo 
can not be explained entirely by M87. In the following we keep the aph = 2.3 results, i.e. M87 contribution of 4.5+0.7 
10-2 c s~\ and the thermal model prediction of 0.3+0.210-^ c s^^ to the PDS 20 - 80 keV band. 



